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PASSIVE DAMPING CONCEPTS FOR SLENDER COLUMNS IN SPACE STRUCTURES 

By 

Zla Razzaq 1 and Rajendra K. Ekhellkar 2 
ABSTRACT 

An experimental and theoretical study of three different passive damp- 
ing concepts Is conducted for a slender member with partial rotational end 
restraints. Over a hundred full-scale natural vibration experiments were 
conducted to evaluate the effectiveness of mass-string, polyethylene tubing, 
and chain damping concepts. The damping properties obtained from the exper- 
iments were used in the approximate analyses based on the partial differen- 
tial equation of motion for the problem. The comparison of the experimental 
and the theoretical deflection-time relations shows that the velocity-depen- 
dent damping model used In the theory Is adequate. From the experimental 
results, the effect of end connection friction and induced axial forces on 
damping Is Identified. The definition of an efficiency Index is proposed 
based on the damping ratio and the mass of a given passive damping device. 
Using this definition, the efficiencies of the three damping devices are 
compared. The polethylene tubing concept resulted into a low damping effi- 
ciency. Although the chain damping concept provides a significant amount of 
damping and a rapid decay of member vibration. Its mass-based efficiency 
Inuex Is much lower than that of the string-mass concept. 

1 Associate Professor, Department of Civil Engineering, Old Dominion Univer- 
sity, Norfolk, Virginia 23508. 

z Graduate Research Assistant, Department of Civil Engineering, Old Dominion 
University, Norfolk, Virginia 23508. 
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1. INTRODUCTION 


1.1 Preliminary Remarks 

Members with large slenderness may be used to form large space struc- 
tures for use In NASA's related space missions Tubular members with slen- 
derness ratios of the order of a thousand may be needed due to the optimiza- 
tion constraints imposed by efficient packaging and deployment requirements. 
A desirable feature of the overall response of a space frame would be its 
ability to transfer the global energy very rapidly, due to an imposed per- 
turbation, to the individual members which may then dissipate the energy 
locally. If this can be achieved, the individual members may be expected to 
act as possible "energy sinks" for the system. Thus, a fundamental problem 
which one is faced with is that of the identification of potential damping 
concepts which may be used to absorb passively all or a large portion of the 
member energy. 

Two types of damping approaches are used currently by engineers for 
absorbing vibration of energy in various kinds of structures (4). These are 
classified as "active" or "passive" damping systems. A damping system is 
active if there- is an external energy supply to reduce to an admissible 
level the vibration characteristics at chosen points of a structure, within 
a given band of frequency or an interval of time for a predicted class of 
excitation. Active systems can generate local forces, which are related to 
distant motion variables or external command signals. Passive systems can 
dissipate energy, or temporarily store it and then return it. Passive damp- 
ing elements generate active or dissipating forces, which are related to 
local motion variables. 

Active damping systems are used where the operating conditions vary so 
widely that variations in the strategy of control are frequently called for. 
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Among the advantages of passive damping systems are their simplicity, con- 
tained weight, low cost, and high reliability. 

i 

This research report presents the outcome of an experimental Investi- 
gation of three different passive damping concepts for a s tender member with 
partial rotational end restraints. Theoretical results based on approximate 
analyses are also presented and compared to those obtained experimentally. 
Partial restraints are used since a member In a real space structure seldom 
possesses pinned end or fixed end conditions. The end connections and the 
adjoining members in a space structure provide partial rotation'*! re- 
straints. 

1.2 Literature Review 

A review of che existing literature indi-ates that a rather limited 
amount of research has been conducted in the past on passive damping of 
slender members. 

Reference 7 describes how a support structure for antennas and solar 
panels can also be used for damping mechanical vibrations. Figure 1(a) 
shov.s the damping support consisting of two thin-wall»d aluminum tubes, one 
of which can slide within the other under vibration loads. Silicone grease 
was used as the vibration-damping mediun in between the two tubes. Figure 
1(b) shows a coil spring attached to the outer tube to transmit the static 
loads imposed on the structure. Under vibration, the coil moves with the 
outer tiie thereby employing viscous drag of the silicone grease to effect 
damping. 

Razzaq, Voland, Bush, and Mikulus (1) presented an experimental and 
theoretical study of the natural vibration of a passively damped uniform 
slender tubular member with partial rotational end restraints as shown in 
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Figure 2. A damping concept consisting of a string-mass assembly was ex- 
plored In a few tests, in addition to the structural damplnq present In the 
member. Lead shots were used as masses. Tests were conducted with no lead 
shots; with one or three lead shots at member mldhelqht; and with three 
equidistant lead shots. late-al vibration tests were conducted In the 
absence of any externally applied axial load. The three-lead conflouratlons 
both provided considerably greater damping than the single lead shot 
damper. 

In the experiments reported In Peference 1, the effect of end connec- 
tion friction was not quantified. Also, there was some axial compress'on of 
•j-'known magnitude as Indicated by the experimental frequencies (average * 
.•M Hz' which were generally somewhat lower than the expected theoretical 
frequencies (average ■ 3.58 Hz). Oue to these factors, and possibly the 
attendant three-dimensional effects, the reported damping ratios were hiqher 
than i-nse found In the present Investigation. 

In Reference 2 there Is a brief outline of the outcome of viscoelastic 
passive damping tests conducted on skeletal space structures. Two huae 
trusses ware erected on the ground for vibration testing. Flqure 3 shows a 
60 ft. high truss subjected to 1 Hz vibration. Aluminum tebinq is being 
used as horizontal members of the test structures. The diagonals of one are 
constructed of plexiglass and the other of Lexan. Damping ftedules on the 
diagonals contain viscoelastic material for passive application and result 
In an Increase of 3. 8X to 6X damping. Extensive studies have been conducted 
In the past on the application of viscoelastic damping materials for plate 
ar d shell-like structures (6). 

In this research report, a study of three different passive damping 
concepts for slender tubular members is presented. 


1.3 Problem Statement 


Figure 2 shows schematically a slender member of length L with a 
hollow circular cross section. The outer diameter of the section Is Do, 
the wall thickness Is t 0 , and the material of the member Is mild steel. 
The sceel member Is used since very slender graphite composite tubes which 
may be used In space structures are not yet available. The member ends are 
provided with connections which possess partial rotational restraint charac- 
teristics. No axial or lateral movement of the member ends Is permitted. 

The problem Is to identify potential passive damping concepts to absorb 
the energy of natural flexural vlbratluii of the member, and to study the 
comparative effectiveness of each concept. The vibration Is Induced by 
plucking and releasing the member at Its midheight. 

1.4 Objective and Scope 

The following are the main objectives of this study. 

1. Identification of potential passive damping concpets for slender 
structural members. Specifically, the following damping concepts 
are Investigated: 

a. Mass-String 

b. Polyethylene Tubing 

c. Chain 

2. Evaluate the suitability of theoretical analyses by comparison to 
the experimental results. 

3. Evaluate the damping efficiencies of the various damping devices. 
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1.5 Assumptions and Conditions 

The following assumptions and conditions have been adopted In this 
study: 

1. The deflections are small. 

2. The material of the member Is linearly elastic. 

3. Only planar vibration of the member Is considered. 

4. Damping force Is opposite but proportional to the velocity at any 
location along the length of the member. 

5. The member Is tested under 1-g conditions at normal room temper- 
ature. 

6. The effect of secondary Induced forces such as the variable axial 
tension or compression developed In the member during vibration Is 
considered to be small and, therefore, neglected In the theoretical 
analysis. 
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2. THEORETICAL BACKGROUND 


2.1 Governing Equations 

The differential equation of equilibrium governing natural vibration of 
a flexural member (Reference 9) can be written as: 

ei£a ♦ p iis ♦ cii ■ 0 (1) 

3x 4 3t 2 3t 

In which 

E ■ Young's modulus of the member material, 

I ■ Moment of Inertia of the member section, 

p ■ Mass per unit length of the member, 

C ■ Damping coefficient, 

w(x,t) ■ Lateral member displacement. 

Equation 1 may be used for the analysis of the member shown In Figure 2 as 
long as the Induced flexural vibrations are small and asstmlng that the 
passive damping devices will provide viscous damp.ng. The validity of the 
latter assumption is Investigated using the experimental data and explained 
in Chapter 3. 

The boundary and initial conditions for the problem in Figure 2 are 
given In Reference 1 and are sunmarlzed here: 



w(o,t) « 0 

(2) 


w(L,t) a 0 

(3) 

EI !!i 

(o.t) ■ K ll (o,t) 

(4) 

3x 2 

3 X 
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El (L.t) - -K — (L.t) 
3>? 3x 


(5) 


3 w 
3t 


(x,o) ■ 0 


( 6 ) 


w(x,o) - 0 (K,EI,L,x) (7) 

Equation 2 to 5 represent the boundary conditions whereas Equation 6 
and 7 are the Initial conditions. Equation 7 expresses the condition that 
at zero time, the member deflected shape is dependent upon K, El, L, and x. 

2.2 Semi-Analytic Solution 

Based on the separation of variables approach, the solution to Equation 
1 may be expressed as: 


r 


i 


w(x,t) » W(x) T(t) 


( 8 ) 


Substituting this equation into Equation 1 results in the following two 
ordinary differential equations: 


d^W 

dx 4 


- XW = 0 


( 9 ) 


d Z T A C dT x , r n 
+ _ + ur T =» 0 

dt 2 p dt 


( 10 ) 


in which x is a parameter to be determined, and u is the undamped clr- 


( 


i 


e* v ' 


. ■*,+ ~ •*- •** **»•?•* ' ' 


*1 


:u1ar frequency given by: 


/e IX 
u ■ / — ■— 

V p 


( 11 ) 


In Reference 1 the following expression for W(x) Is adopted: 


W(x) - Ai (sin — + JL. (i . cos — )] 
L 4»EI L 


2*xi 


( 12 ) 


If the mid-height displacement of the member Is A 0 , then Equation 12 
gives: 


A i * 


1 + 


KL 

2*EI 


(13) 


Invoking Galerkln's orthogonality criterion, the following expression for X 
is obtained (Reference 1): 


X . ( n )' < 12(*EI) 2 + 80EIKL + 1 2 ( KL ) 2 

(L)“ 12(*EI) 2 + 32EIKL + 2.25 (KL) 2 


The critical damping factor Is given by (Reference 9): 


c c - 2 VTlpX 


(14) 


(15) 


The damping coefficient C in Equation 10 is given by: 


8 




(16) 


In which c, the damping ratio, Is obtained using the logarithmic decrenent 
method (5). The closed-form solution of Equation 10 for an underdamped 
system is (9): 


T(t) - e " ct/2p (Fi cos 8t + F 2 sin o t] (17) 

where Fi and F 2 are constants of integration and 8 is given by: 

8 « 0.5 yj 4w 2 - (c/p) 2 (18) 

Using Equations 6 and 8 one obtains: 

. o (i« 

dt 

which gives 

F 2 = (20) 

P 8 

when substituted into Equation 16. From Equations 12, 17 and 20 one ob- 
tains: 

_ ct 

w(x,t) = A [ sin Hi + (l - cos ^Hi)] e 2p [cos 8t 

C L 4tr E I L 


( 21 ) 


£ 

+ — sin 8t] 
pB 
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l In which A Is a constant, and Is found to be equal to Ai since the 
t * 

; Initial deflection at x » 0.5L Is a 0 . Equation 21 can be used for a 
| direct computation of the member lateral deflections as a function of time 
i and the various Input parameters. 

3.3 Finite-Difference Solution 

Using central finite-difference expressions with second order trunca- 
tion errors (Reference 8), Equation 1 can be written as follows (Reference 
1 ): 


^ ^ j ^ ^ j ) + 6w(x^,tj) 


4w( x i+i» tj) + t^] + “ — rr [w(x <t t 1+1 ) 


'i+2* ‘J 


(At) 




2w(x^, t j ) + w(x.j, t^_j)] + [ w ( x ^» tj + j ) 


- w ( x r tj)] = 0 


( 22 ) 


The above equation can be abbreviated as follows: 


E]_ 

h- 


' 4 Vi. j * 6 "i, j - 4 Vi, j 

* "i+2, j* * ^ W 1, j+1 ' 2u i,j * 
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in which 


* ^"un - w i,j ! ;° 


(23) 


h 

At 

x i 

m 

j 


max 


■ panel length along the x axis of the member, 

■ small time interval, 

■ ih for each i » 1,2,3.... M, 

■ j(At) for each j 0,1, 2, 3,..., j max , 

■ number of equidistant nodes over [0,1], 

« maximum desired j values. 


Similarly, the boundary conditions 2 to 5 in the difference form becomes: 



w . » 0 
o,J 

(24) 


W_ 2 ■ 0 

m, j 

(25) 

EI_ + K 
h 2 

El K _ 

w -i,j 7 w u 3 0 

(26) 

IL*!L 

h 2 

W M-l,j + J W M-l,j * 0 

(27) 


Applying Equation 23 at i ® 1,2,3 ... (M-l) and invoking Equations 24 to 27 
leads to the following matrix equation: 

{ w i , j+i) ■ Nij) + C2 ^ w i,j-i^ (28) 


in which 


11 


b3 + 


C2 M b 3 C 1 


b, ■ — 


n _ k 

h ’ 2 


EI_ + K 
h 2 


b - • — 


where [ r ] Is the following symmetrical square matrix of orde 


Y 2 f 3 


Y3 Y 3 


Y2 Y3 


[r] = 


Symmetric 


Y4 


.1 




j In which 


Yi ■ 6bj - 2b3 - bi* - bj b s (36) 

Y 2 ■ 6b x - 2b 3 > h, (37) 


Y3 * -4bi 


(38) 


Y4 * bj (39) 

Equation 28 can be used to calculate the deflection if (w^ and 

(w. . .) are known. To avoid a negative time interval owing to the use of 

1 * I t J 

central finite-difference scheme when starting the recursion indicated in 
this equation, the following special forward start-up difference equation 
for Wj ^ is derived in Reference 1: 

w 1,l = bsl,, i-2 “ 4t)6 *i-l + + 1) + b6 *1+2 ( 40 ) 


in which 


b6 . .Miji. 

2h 4 p 

The initial static deflected shape in Equation 7 may be determined 
using any one of the several classical structural analysis, an approximate 
deflected shape (Reference 1) for the member due to a midheight displacement 
A 0 at t=0 is taken in the following form: 
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(41) 


W(x) Is given by Equation 12. 

The following steps sunmarlze the procedure for determining the approx- 
imate deflection time response of the member: 

1. Define input data: L, K, p, E, I, M, AT, and j . 

2. Calculate the constants b t , b 3 , b 4 , b 5 , q , c 2 , using Equations 29- 
34; Yj to y 4 using Equations 36-39; and formualte [r] using Equa- 
tion 35. 

3. Calculate X using Equation 14. 

4. Calculate the critical damping facior C cr using Equation 15, and 
the damping coefficient C using Equation 16. 

5. Set j = 0 and t » 0. 

6. Specify the midspan deflection A 0 and generate the initial de- 
flection vector (w. n ) using Equation 41. 

7. Set j ■ 1, t = At, jnd calculate {w^ ,} using Equation 40. 

* 

8. Increment j by 1. 

9. Set t = j (At) . 

10. Calculate (w^ using Equation 28. 

11. If j < j max ! go to 9. 

12. Stop. 

In the above procedure, At was taken as 0.001 sec. based on a prelim- 
inary time-convergence study. Also, dividing the member into 11 equal 
panels (M=12) was found to be appropriate for the nunerical study present- 


14 


ed In the following section, as well as, for comparisons to the experimental 
data given In Chapter 3. 

2.4 Numerical Results 

For the deflection-time analysis of the member shown in Figure 2, com- 
puter programs were written based on the procedures given in Section 2.2 and 
2.3, respectively. In this section, a numerical study based on these 
computer programs is presented with the following input data: L 3 12 ft.; 
D 0 * 0.5 In.; t 0 = 0 065 in.; K 3 2,230 lb-in/rad; p = 651 x 10 -7 lb-sec 2 / 
In 2 ; E 3 30 x 10^ psi; and A 0 = 0.366 in. This data corresponds to the 
vibration tests described in Chapter 3. 

Figure 4 shows deflection-time plots obtained using the semi-analytic 
and the finite-difference methods with the above mentioned input data and 
C 3 0.006. The undamped natural frequency using the finite-difference ap- 
proach was observed to be 3.53 Hz whereas that obtained from the semi- 
analytic approach ’./as 3.73 Hz. Table 1 shows a comparison of the deflec- 

tion-time envelope ordinates obtained from Figure 4 based on the two methods 
for t 3 1 to 6 seconds. A comparison of both positive and negative enve- 
lope ordinates is made. In this table, App and a^ represent the ordi- 
nates of the envelopes based on the finite-difference, and the semi-analytic 
procedures, respectively. The ratios App/A^ are also given in the table. 
As evident from the table, the ratio. App/A^ ft increases as time increases. 

The maximum deviation between the two methods is found to be about 2.50 
percent, at t = 6 seconds. 

Figures 5 and 6 show the envelopes for various values of c in the 
range from 0.0 to 1.0 and obtained by the finite-difference and the semi- 
analytic procedures, respectively. A comparison of the curves in Figures 5 


i5 



and 6 shows that the maximum different between the results of the two meth- 
ods Is of the same order of magnitude :.s discussed above for Figure 4. In 
addition. It was observed that this difference diminished gradually with an 
Increase In c. Computer Programs and their sample outputs are given In 
Appendix A. 
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3. tXPERIMENTAl STUDY 


I 

i 

3.1 Specimen and End Fixtures 

The experimental study consisted of conducting natural vibration tests 
on a tuhular steel member. The tests Mere performed both with and without 

i 

: passive damping devices. The member ends were provided with specially de- 
i signed base plates which In turn were bolted to external fixed supports. 

The connections between the member ends and the base plates provided partial 
■ rotational restraint. The tubular member was 12 ft. long with Dq ■ 0.5 In, 
j and tg ■ 0.065 In. 

Figure 2 shows schematically the member tested. The member end con- 
nection details are shown In Figure 7. Each end connection was a two piece 
assembly. One piece was a steel blade 1-1/2 in. long with a 3/16 x 1/2 In. 
cross section welded to a steel base plate. The outer piece was a steel 
clevis designed to fit snugly overall but 1/4 in. of the blade and fastened 
to It by screws fitted through two drilled holes. The exposed 1/4 In. of 
the blade, thin In one direction (3/16 In.) but thick In the other (1/2 in.) 
provided partial rotational restraint in one plane of vibration and consid- 
erably greater rotational restraint in the orthogonal plane. This design 
provided the member a preferred direction of vibration, that is, about the 
weaker axis of the blade. The top end of the clevis was a 1 In. threaded 
rod which screwed by a distance S * 5/8 In. Into the specimen. Since the 
primary purpose of the experimental Investigation was to study the effect on 
C of each passive damping device under consideration, the connection 
friction was reduced considerably by means of a "jam nut" screwed on to the 
threaded clevis rod and butted against each member end, as shown in Figure 
7. The jam nuts were used for all of the tests reported in this section 
except for an initial set of tests which were conducted for the purpose of 
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comparison. The upper base plate was bolted to a heavy bracket which In 
turn was bolted to a backstop, and the lower base plate was bolted to the 
floor directly beneath the upper one. 

3.2 Specimen Initial Imperfection 

The tubular member used In the experimental part of this study had an 
Initial Imperfection. The member deviation at right angles to the plumb 
line was measured using a theodolite placed at a distance of 7.2 ft. from 
the member. Table 2 summarizes the specimen initial Imperfections measured 
at 1 ft. Interval:. The second and the third colunns In this table present 
the measured deviation of the member centerline relative to the plumb line 
in the plane of vibration and In the orthogonal plane. As seen from the 
deviation at 12 ft. from the top of the member, the member wa: out-of-plumb 
by an absolute amount of 0.10 in. In each plane. The fourth and fifth col- 
umns of Table 2 represent the actual Imperfections in the two planes calcu- 
lated from the deviation measurements. Figures 11(a) and 11(b) show graphi- 
cally the data give In Table 2. In these figures, the curved solid lines 
from A to B represent the actual member centerline. The plumb lines are 
represented by the chain lines BA' and BA", and the straight dashed lines 
from A to B are used to define the actual imperfection at 1 ft. Intervals. 
For proper visualization of the spatial imperfection, the probe location Is 

also identified In these figures. The maximum imperfection values In the 
two planes were found to be 0.070 In. and 0.378 in., respectively. Such 
small imperfection is not expected to affect the natural frequency of the 
member. 
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3.3 Test Procedure and Instrumentation 

The end restraint rotational stiffness K about the minor axis was eval- 
uated experimentally for both top and the bottom end connections using the 
procedure given In Reference 1. Twelve experiments were conducted on each 
end connection and the moment-rotation relations obtained. Figure 8 shows a 
typical moment-rotation relationship for the bottom end connection. All 
such relationships were found to be linear. The average value of K was 
found to be 2,230 lb-ln/rad. Similar tests were conducted to obtain major 
axis end restraint rotational stiffness. Figure 9 shows a typical major 
axis moment- rot at Ion relationship for the bottom end connection. The aver- 
age major axis rotational connection stiffness was found to be 7,280 lb- 
ln/rad. Since the slender member under Investigation Is vibrated about the 
major axis, the applicable K valu? for the results presented later In this 
research report Is 2,230 lb-in/rao. 

The Instrumentation used In the test consisted of a proximity probe, 
vibration Instrunentatlon, and a deflection-time curve plotter. Figure 10 
shows the proximity probe mounted on a bracket at member midheight in the 
plane of vibration. To Induce natural vibration, the specimen was pulled 
away from the proximity probe by 0.366 In. with the help of a nylon string 
and then released by burning the string with a match stick. 

The midheight deflection versus time curves were obtained automatically 
by means of an on-line plotter. Damping ratios were computed from each 
deflection-time plot using the logarithmic decrement method (9). Each test 
was repeated three time to obtain a good average value for the damping 
ratio. The details of the Instrumentation and the initial setup are ex- 
plained in Appendix B. 


19 



3.4 String-Mass Damping Concept 

In Reference 1, vibration tests were conducted on an Imperfect tubular 
steel member with partial rotational end restraints. A passive damping 
concept consisting of a string-mass assembly was explored In a few of the 
tests, In addition to the structural damping present In the system. Lead 
shots were suspended Inside the tubular member by means of a nylon string. 
In addition to the baseline tests with no lead shots, experiments were con- 
ducted with one and three lead shots suspended inside tne member. The tests 
were conducted without any jam nuts at the member ends. The test results 
reported In Reference 1, therefore, represented the combined effects of the 
connection friction, the string-mass dampers, inherent structural damping, 
and other Influences on the Indicated damping characteristics of the member. 
The other Influences Included the presence of the Initial imperfection, 
small Induced axial forces, and three-dimensional movement of the member due 
to a spiral-like motion emanating from the threaded end connections In the 
absence of jam nuts. In the experiments reported herein, the installation 
of jam nuts was seen to result In a considerable reduction of the three- 
dimensional member movement. 

As a part of the present research, the string-mass concept was studied 
extensively both with and without jam nuts to identify the effects of the 
connection friction and the number of lead shots on the member damping. 

Besides the baseline experiments on the member without any damping de- 
vices, tests were conducted with a total of twelve different string-mass 
configurations shown in Figures 12 and 13. Figures 12(a) through 12(h) show 
1, 2, 3, 5, 7, 9, 15 and 21 lead shots suspended inside the member by a 
nylon string. The lead shots were attached to a nylon string with a 0.2 cm. 
clear spacing between them. The gap between the lead shots allowed flexible 
movement of the configurations during vibration. Figures 13(a) through 
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13(c) show 2, 3, and 21 lead shots suspended equidistantly Inside the member 
with spaclngs of 1/3, t/4, and 1/22, respectively. The lead shots were 
mounted symmetrically about the member midheight and the nylon string did 
not extend beyond the lower most lead shot. Each lead shot has an average 
diameter of 0.254 cm. and weighed 0.47 gm. The nylon string weighed 0.01417 
gms. per ft. As shown In Figure 14, the nylon string was tied to an 
Internal screw which In turn was screwed Into the top threaded clevis 
rod. 

3.4.1 Performance of String-Mass Concept with Lead Shots at Midheight 

A series of natural vibration tests were conducted on the member shown 
In Figure 2 and provided with the string-mass configurations shown In Figure 
12 both In the absence and presence of jam nuts. In addition to the base- 
line natural vibration tests In the absence of passive dampers, the experi- 
mental results for the eight string-mass configurations shown In Figure 12 
without jam nuts are summarized in Table 3. Each configuration was tested 
three times to obtain a good average. Tests 1-3, 4-6, 7-9, ... were per- 
formed with 0, 1, 2, ... nunber of lead shots, respectively, as Indicated In 
the first two columns of Table 3. The experimental plots of the midheight 
deflection. A, versus time, t, for Tests 1, 4, 7, 10, 13, 16, 19, 22, 
and 25 are shown in Figures 15 through 23, respectively. These are repre- 
sentative plots for Tests 1-27. The values of the experimental natural 
frequency, f , and the damping ratio, c, as well as their averages for 
these tests are given In columns 3-6 of Table 3 and were computed using the 
deflection-time plots. The average f values for the various tests ranged 
from 3.84 to 3.98 Hz. The average c values Increased substantially with 
an Increase In the number of lead shots and ranged from 0.00243 to 0.00693. 

Table 4 summarizes the experimental results for Tests 28-54 in the 
presence of jam nuts and correspond, respectively, to Tests 1-27 as far as 
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their string-mass configurations are concerned. Representative A-t plots 
for the tests referred to In Table 4 are given In Figures 24-32. The aver- 
age f fi values for the various tests ranged from 3.74 to 3.91 Hz. The 
average c values ranged from 0.00103 to 0.00615 with the upper limit cor- 
responding to 21 lead shots. 

In the experiments with 21 lead shots, the value of c was found to 
vary slightly with time. This variation diminished as a decreasing number 
of lead shots was used. 

3.4.2 Effect of Connection Friction on Damping Ratios 

Figure 33 shows the relationships between the damping ratio, c, and 
the number of lead shots, n, both with and without jam nuts. The two 
curves shown In this figure are plotted using the data given In the second 
and the sixth columns of each Tables 3 and 4. The solid curve for the tests 
with jam nuts and the dashed curve for those without are fairly parallel to 
each other Indicating a nearly constant effect on c of the connection 
friction and other possibly less dominant Influences described earlier in 
this research report. The average constant difference between the c 
values of the two curves is 0.0012. Also, the shape of the c-n curves in 
Figure 33 shows a rapid initial Increase in ;. The rate of increase in c 
decreases gradually with an Increase in n. 

3.4.3 Performance of String-Mass Concept with Equidistant Lead Shots 

In order to study the Influence of the location of lead shots on damp- 
ing, a series of vibration tests were conducted on the member provided with 
the string-mass configurations with equidistant lead shots in Figure 13. 
The end connection jam nuts were used In all the tests. 

Table 5 summarizes the experimental results for Test 55-63 for 2, 3, 
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and 21 lead shots spaced equldlstantly. Representative A-t plots for 
these tests are given In Figures 34-36. The respective average f values 
were found to be 3.93, 3.89, and 3.98 Hz. The corresponding average c 
values were 0.00228, 0.00258, and 0.00462. 

3.4.4 Effect of Location of Lead Shots on Damping Ratio 

A comparison of the damping ratios for the 2 and 3 lead shot tests 
given In Table 4 to those In Table 5 Indicate that they are practically the 
same. This conclusion, however, does not hold for the 21 lead shot tests as 
evidenced by the significant reduction of c from 0.00615 for the midheight 
location to 0.00462 for the equidistant case. Figure 37 shows a comparison 
of the corresponding A-t riots. From these results, it Is concluded that 
the midheight location for the lead shots provides better passive damping 
than the lead shots spaced equldlstantly. 

3.5 Effect of Induced Axial Force on Member Frequency and Damping Ratio 

An axial tensile or compressive force can be Induced In the member by 
turning It about the longitudinal axis (x axis In Figure 2) In the clockwise 
or counterclockwise directions, respectively. A series of member vibration 
tests were conducted with Induced axial forces and in the absence of any 
damping device. Jam nuts were used In all the tests. Table 6 summarizes 
the results for Tests 64-75 for four different values of axial force, P. 
The values of P recorded in this table were not measured experimentally 
but were determined indirectly by solving Equation 28 of Reference for P 
in terms of the natural frequency f, and other parameters as follows: 

P - L (u 2 - 4 f2 L2 p U3 ) (42) 
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n which 


Wl - 12(itEI) 2 + 32EIKL + 3(KL) 2 

(43) 

, (12(*EI) 2 + 80 EIKL + 12(KL) 2 ] 

(44) 

l 2 


, • 12(uEI) 2 + 32 EIKL + 9/4 (KL) 2 

(45) 


The undamped natural frequency for each test can be obtained using the fol- 
lowing expression (Reference 9): 




(46) 



The average fe and c values given In the last two columns of Table 6 
: W ere used to find f through Equation 54. The P values were then com- 
puted using Equation 42 and entered In the second column of Table 6. The 
negative and the positive signs for P In this table are used for tensile 
and compressive forces, respectively. For Tests 64-66, the member was ro- 
tated In a clockwise direction until It could no longer be rotated further 
manually. The computed P was found to be a tensile force of -155.0 lb. 
Similarly, other amounts of member rotations resulted in axial forces of 
-19.8, -4.5, and +10.8 lb., respectively, for Tests 67-69, 70-72, and 73-75, 
as indicated in Table 6. 

As expected the average member frequency f £ Increases from 3.83 Hz to 

i 
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1 6.64 Hz as the tensile force Is Increased from -4.5 1b. to -155.0 1b., 
whereas damping ratio c experiences some reduction, as Indicated In Table 
6. Similarly, the average f fi decreases, attaining a value of 3.42 Hz when 
an axial compressive force of +10.8 lb. Is Induced, with a significant In- 
crease In c, as given In Table 6. The representative a-t plots for the 
tests are given In Figures 38-41. 

Tests 70-72 resulting In P » -4.5 1b. deserve special attention. In 
these tests, the member was adjusted by axial rotations with a view to ob- 
tain a zero axial force. This was attempted using an experimental trial and 
error procedure as follows. When a small axial compressive force was In- 
duced, the member appeared slightly bent by visual observations. On the 
other hand, when a small axial tensile force was Induced, the member fre- 
quency showed a corresponding slight Increase measured from a a-t plot. 
For a "zero" point, that Is, for no expected axial force, the approximate 
position of the member between the compressive and the tensile states was 
determined by trial and error. The vibration tests were then calculated 
with the assunptlon that the axial force would be nearly zero. The computed 
axial force as presented In Table 6, however, turned out to be -4.5 lb. The 
presence of this force may be attributed largely to the approximate nature 
of the trial and error procedure as well as the unavoidable tension induced 
In the member during vibration since the member ends are not allowed to move 
axially. From these results. It Is concluded that such small amounts of 
unavoidable axial forces were present in all of the tests conducted in this 
research. 
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3.6 Polyethylene Tubing Damping Concept 

A series of natural vibration tests were conducted on the member shown 
in Figure 2 and provided with various types of polyethylene tubing over Its 
entire length. All the tests were conducted with jam nuts. The following 
three tubing arrangements were Inserted in the member: 

(a) A low stiffness polyethylene tube with an outer diameter of 0.33 
In., a wall thickness of 0.04 In., and weighing 0.0254 Ib/ft. 

(b) A stiff polyethylene tube with an outer diameter of 0.33 In., a 
wall thickness of 0.04 In., and weighing 0.0140 Ib/ft. 

(c) A low stiffness polyethylene tube with an outer diameter of 0.25 
In., a wall thickness of 0.04 in., and weighing 0.00875 lb/ft Inserted In 
the stiff polyethylene tube (b) and the resulting assembly Inserted Into the 
member. 

Table 7 sunmarizes the experimental results for Tests 76-81 based on 
the above mentioned polyethylene tubing concepts. The deflection time plots 
for Tests 76, 79, and 82 are shown In Figures 42-44, respectively. The 
maximum c was found to be 0.0022 for the tubing arrangement (a) and is 
nearly the same as for Tests 34-36 (Table 4) which were conducted with two 
lead shots at member midheight. Clearly, the polyethylene tubing concept is 
far less effective than the string-mass concept in providing passive damp- 
ing. 

3.7 Chain Damping Concept 

A "bright zinc chain" was suspended inside the tubular member by a 
nylon string to absorb the vibration energy through collisions as well as 
frictional forces between the chain links. Figure 45 shows the photograph 
of a portion of the chain used for the experiments. The chain weighed 0.05 
Ib/ft and its major lateral dimension was approximately equal to the inner 
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diameter of the tubular member. Chain*; of 1, 2, 3, 5, 7, 9, and 10 ft. 
lengths were tested. Each chain was situated symmetrical 1y about the member 
midheight. Natural vibration tests were conducted In the presence of end 
jam nuts. 

Table 8 summarizes the experimental frequencies for Tests 85-105 for 
chain lengths ranging from 1 to 10 ft. Deflect Ion- time plots for Tests 85, 
88, 91, 94, 97, 100, and 103 are given In Figures 46-52, respectively. A 
noteworthy characteristic of the chain damping approach Is that the damping 
ratio c obtained from It Increases considerably with time as the chain 
length is Increased, unlike the mass-string or polyethylene tube concepts In 
which c remains constant. Table 9 presents the c values for Tests 85- 
105 at average times of 0 to 6 seconds with one second intervals. The cor- 
responding average c values for various chain lengths are given In Table 
10. The data given In Table 10 is also shown graphically In Figure 53 Indi- 
cating the time dependence of c In the form of curves relating c to the 
average time, t . From this figure, it Is clear that c Increases rapidly 

a 

both with an Increase In the chain length as well as the time of vibration. 
The best results are obtained for the 10 ft. long chain for which c In- 
creases almost exponentially with time. The average damping ratio, 

av 

versus the chain length relationship Is shown In Figure 54. The C av 
values used in this plot were obtained by taking the averages of z values 

/ 

/ given In Table 10 for the average times from 0 to 6 seconds for each chain 

length. 
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4. COMPARISON OF RESULTS AND DAMPING EFFICIENCIES 



1 Theory Versus Experiment 

Two different theoretical solutions, namely, the semi -analytic and the 
inite-difference procedures, for predicting the natural vioration response 
F the member shown in Figure 2 were outlined In Chapter 2. Experimental 
esults for the string-mass, polyethylene tubing, and chain damping concepts 
'or a total of 105 tests were presented in Chapter 3. In this section, a 
omparlson of the experimental results is made to the theoretical ones for 
'ests 52, 76, and 103 corresponding to the maximum c values obtained for 
:he three concepts, respectively, which were conducted in the presence of 
the jam nuts. 

Figures 55-57 show the comparison of the experimental and the semi- 
analytic deflection-time (A-t) plots for Tests 52, 76, and 103, respec- 
ti«e1y. for the three concepts. The natural frequency from the semi-ana- 
lytic procedure was found to be 3.73 Hz. Figure 55 shows the A-t plots 
for Test 52 corresponding to 21 lead shots at the member midheight. Two C 
values equal to 0.0085 and 0.00615 obtained from the experimental A-t plot 
for the first, and all of the suDsequent cycles, respectively, were fed into 
the semi-analytic procedure. A higher first c value is needed for a prop- 
er accommodation of the initial nonlinear part of the A-t relationship. 
The experimental natural frequency was found to be 3.78 Hz. Similarly, 
Figure 56 ror Test 76 shows the comparison of the A-t plots for the poly- 
ethylene tubing concept (a), with c « 0.0096, and 0.00220, respectively, 
for the first, and all of the subsequent cycles. The experimental natural 
frequency was found to be 3.96 Hz. Figure 57 for Test 103 shows the com- 
parison of the A-t plots for the chain damping concept. The experimental 
frequency was found to be 4.01 Hz. As mentioned in Chapter 3, the c 



si ues for the chain damping concept varied considerably with time. Thus, 
ir a proper comparison of the experimental and the theoretical A-t plots, 
ie following variable parabolic c expression was developed using the 
/erage c values at 0, 3, and 6 seconds, for the 10 ft. chain as given In 
able 10. 

4 - 0.001289 t 2 - 0,002270 t + 0.008570 (47) 

Jtrictly speaking. Equation 21 In Section II Is applicable only when the 
damping coefficient C, related directly to c via Equation 16, Is a con- 
stant. Obviously, Equation 47 leads to a time dependent C expression. 
However, Equation 21 was used directly with the C expression based on 
Equation 47 to obtain the theoretical A-t curve shown In Figure 57 and was 
found to be in fairly good agreement with the experimental curve. This may 
be attributed principally to the fact that the variation In C with t is 
not really drastic at least during the first 3 seconds; it varies from 
0.00852 to 0.01320. The subsequent variation in 5 is more rapid, from 
0.01320 to 0.04090 in 3 additional seconds, however, it did not also result 
In to a very significant deviation of the theoretical A-t plot relative to 
the experimental one except In the region of the fairly low-amplitude vibra- 
tion toward the tail end of the A-t relationship. 

Figures 58-60 show the comparison of the experimental and the finite- 
difference A-t plots for Tests 52, 76, and 103, respectively, for the 
three concepts. The natural frequency from the finite-difference procedure 
was found to be 3.53 Hz. As seen from these figures, the experimental and 
the finite-difference A-t plots are in good agreement. The same c values 
were used for the finite-difference curves as were used in the semi-analytic 





ocedure, Including Equation 47 for chain damping. It should be pointed 
it that a time-dependent c can readily and accurately be accounted for In 
ie finite-difference solution owing to its discretized nature relative to 


.2 Efficiency of Various Damping Concepts 

The main purpose of a passive damping device is to damp out the deflec- 
lons of a member In minimum possible time. A damping device would be 
onsldered efficient If It would achieve this goal with a minimum of Its own 
ass. Any kind of a measure of the efficiency of a damping concept, there- 
ore, should Involve a damping property and the mass of the passive damping 
evice. For evaluating the relative merit of the passive damping concepts, 
he following definition of an "efficiency index" is proposed herein: 

n » L :JS. (48) 

m d 

in which c and m^ arc the damping ratio, and the mass of the passive 
damping device, respectively, and c 0 is the damping ratio in the absence 
of any passive damping de-'ice. The subtraction of from c eliminates 
the effect of the inherent structural damping so that only the effect of 
damping due to the damping device itself is reflected in the efficiency 
index. 

Table 11 presents the computed efficiencies of the various damping 
concepts considered in this investigation. The test numbers, the details of 
the damping devices, the weight w d of the passive damping device, the 
average z, and the efficiency index n based on Equation 48 are given in 
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this table. The mass of the passive damping device Is computed using the 
following expression: 

U. 

m d - J. (49) 

9 

In which g * 32.2 ft/rec* 1 . The baseline Cq value Is taken as 0.00243 as 
obtained from Tests 1-3 for computing n values for the string-mass concept 
without jam nuts. For the remaining tests for all three concepts In the 
presence cf jam nuts, Cc Is taken as 0.00103 as obtained from Tests 28-30. 
Also, for the chain damping concept, the average c values given in Table 
11 represent the average of the average damping ratios over 0-6 seconds 
range for each chain length given In Table 10. 

A glance at the n values In Table 11 shows that It varies from 0.12 
tor the tubing concepts (a) and (c) to 19.44 for the string-mass concept 
with two equidistant lead shots. All other efficiency Index values for the 
three concepts lie between these two extremes. For the string-mass concept 
with lead shots at midheight and without jam nuts, n varies between 6.67 
and 17.11, while that with jam nuts varies between 7.59 and 19.28. For 
equidistant - lead shots, n * 5.31 to 19.44. For the polyethylene tubing 
concept, n * 0.12 to 0.20. For the chain damping concept n r 0.56 to 
1.21. Figure 61 shows the experimental relationships between n and the 
number of lead shots for the string-mass concept, and the chain length for 
the chain damping concept. The lower horizontal axis of this figure gives 
the chain length while the upper one defines the number of lead shots, n. 
The upper two relationships are for the string-mass concept with and without 
jam nuts. The three data points for the equidistant lead shots are also 
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shown In this figure. The relationship for the chain damping concept Is 
shown below the string-mass concept curves. A, evident from Figure 61, the 
maximum n value Is obtained for two equidistant lead shots. Also, two 
lead shots give the peak points on each of the n-n relationships ' *'th 
with and without jam nuts. As the number of lead shots Is Increased beyond 
two, the ii value decreases almost continuously. Although the chain 
damping 4-t relat onshlps show a rapid decay of member vibration (see 
Figure 60) as compared to, for example, the 21 lead shot case (Figure 58), 
the mass-dependent definition of the efficiency Index used herein makes the 
chain damping concept appear less attractive. Another Interesting point of 
observation is that the presence of significant end connection friction 
generally reduces the efficiency Index of the string-mass damping concept. 
This can be seen clearly from Figure 61 In which tne n-n relationship for 
the string-mass concept with jam nuts is generally above that for the one 
without jam nuts. 
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S. CONCLUSIONS ANO FUTURE RESEARCH 

5.1 Conclusions 

The following principal conclusions are drawn from this research on a 
partially restrained slender tubular member subjected to natural vibration: 

1. The semi-analytic and the finite-difference solutions are in rea- 
sonable agreement with the experimental results. 

2. The effect of end connection friction on damping ratio is nearly 
constant for various masses of the passive damping devices. 

3. The presence of an axial compressive force In the member increases 
the apparent damping ratio while a tensile force decreases It. 

4. For the string-mass concept, the damping ratio increases gradually 
at a decreasing rate as the number of lead shots is Increased. 

5. Equidistant lead shots result in lower damping ratios as compared 
to the lead shots at member midheight except for the two lead shot 
case. 

6. The polyethylene tubing concept does not provide any significant 
damping. 

7. For the chain damping concept, the damping ratio increases very 
rapidly with an increase in the chain length as well as with time. 

8. Although the chain damping concept provides a significant amount of 
damping and a rapid decay of member vibration, its mass-based effi- 
ciency index is much lower than that of the string-mass concept. 

9. An efficiency Index for measuring the damping performance of the 
various passive damping devices is proposed and shows that the 
string-mass concept is more efficient than the chain damping and 
the polyethylene damping concepts. 
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5.2 Future and On-Going Research 

In the present Investigation, three different types of passive damping 
concepts were studied. Other passive damping concepts should be explored 
with a view to maximize the damping efficiency index as proposed herein. 
These may include electrostatic damping, suspended chambers with oil and 
discs, and other hybrid concepts. Experimental work on some of these con- 
cepts Is In progress. 
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NOMENCLATURE 


Damping Coefficient 
Critical damping coefficient 
Young' r modulus 
Undamped natural frequency 


Damped natural frequency 

Acceleration due to gravity 

Panel lengtn along tne length of the memOer 

moment of Inertia of the memher section 

Ortt a tional stiffness 


Length of the member 

Hass of the passive damping device 

Number of lead shots 

Axial force 

Time 


Weight of the damping device 
Lateral member displacement 

Deflection ordinate based on 


th e finite-difference procedur 


Initial member deflection 
Deflection ordinate based on 


the Semi-analytical procedure 


Small time increment 
Oamping ratio 


Damping ratio without damping device 

Efficiency index 

... ipnath of the member 
Mass per unit teng 
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'able 1. Comparison of Deflection-Time Envelope Ordinates (c ■ 0.006) 



Positive deflection envelope 

Negative deflection envelope 

T ime. .t 
seconds 

4 , 

FD 

Inches 

4 , 

.SA’ 
Inches 

a FD / a SA 

A f0» 

Inches 

a sa» 

Inches 

a FD / a SA 

1 

0.314 

0.314 

0.999 

0.310 

0.311 

0.995 

2 

0.284 

0.287 

0. 988 

0.272 

0.275 

0.986 

3 

0. 243 

0.249 

0.985 

0.242 

0. 246 

0. 983 

4 

0.212 

0.216 

0.981 

0. 205 

0.210 

0.978 

5 

0.181 

0.185 

0.978 

0.176 

0.180 

0.976 

6 

0.156 

0.160 

0.975 

0.154 

0.159 

0.975 
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able 2. Specimen Initial Imperfection 


Distance 
from top 
of member, 
feet 

Deviation, in 

Inches 

Actual Imperfection, In Inches 

In plane 
of vibration 

In orthogonal 
plane 

In plane 
of vibration 

In Orthogonal 
pi ane 

0 

0. 000 

0.000 

0.000 

0.000 

1 

0.000 

0.050 

0.008 

0.108 

2 

0.000 

0.150 

0.023 

0.216 

3 

0. 010 

0.250 

0.075 

0.275 

4 

0. 020 

0.310 

0.070 

0.333 

5 

0.020 

0.420 

0.068 

0.361 

6 

0.020 

0.410 

0.070 

0.360 

7 

0. 010 

0.420 

0.061 

0.378 

8 

0.010 

. 0.400 

0.053 

0.276 

9 

0.000 

0.350 

0.045 

0.250 

10 

-0.060 

0. 300 

0.016 

0.130 

11 

-0. 100 

0.200 

0.008 

0.047 

12 

-0.100 

0.100 

0.000 

0.000 
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ible 3. Experimental Results for String-Mass Concept with Lead Shots 
at Midheight without Jam Nuts 


Test 

No. 

Nutiber of Frequency 

Lead shots f , Hertz 

e 

Damping 
ratio, c 

Average 
f , Hertz 
e 

Average 

C 

1 

3.98 

0.00252 



2 

0 3.95 

0.00256 

3.97 

0.00243 

3 

4.00 

0.00221 



4 

3.95 

0.00273 



5 

1 3.92 

0.00240 

3.93 

0.00278 

6 

3.92 

0.00321 



7 

3.87 

0. 00294 



8 

2 3.98 

0.00344 

3.92 

0.00353 

9 

3.93 

0. 00428 



10 

3.96 

0.00394 



11 

3 3.98 

0.00407 

3.98 

0.00395 

12 

3.99 

0. 00386 



13 

3.92 

0.00439 



14 

5 3.96 

0.00432 

3.93 

0.00429 

15 

3.91 

0.00418 



16 

3.78 

0. 00479 



17 

7 3.91 

0.00461 

3.87 

0.00466 

18 

3.92 

0.00459 



19 

3.92 

0.00487 



20 

9 3.80 

0.00539 

3.84 

0.00515 

21 

3.82 

0.00522 



22 

4.12 

0.00574 



23 

15 3.90 

0.00592 

3.94 

0.00601 

24 

3.80 

0.00637 



25 

3.98 

0.00668 



26 

21 3.82 

0.00683 

3.84 

0.00693 

27 

3.73 

0.00730 






able 4. Experimental Results for String-Mass Concept with Lead Shots 
at Midheight with Jam Nuts 


Test 

No. 

Nunber of 
Lead shots 

Frequency 
f , Hertz 
e 

Danplng 
ratio, ; 

Average 
f , Hertz 
e 

Average 

C 

28 


3.57 

0.00110 



29 

0 

3.85 

0.00100 

3.74 

0.00103 

30 


3.82 

0.00100 



31 


3.83 

0.00140 



32 

1 

3.83 

0.00133 

3.81 

0.00145 

33 


3.79 

0.00165 



34 


3.81 

0.00220 



35 

2 

3.92 

0.00248 

3.74 

0.00227 

36 


3.51 

0.00215 



37 


3.91 

0.00265 



38 

3 

3.89 

0.00245 

3.91 

0.00263 

39 


3.93 

0.00280 



40 


3.81 

0.00275 



41 

5 

3.84 

0.00285 

3.83 

0.00181 

42 


3.84 

0.00285 



43 


3.76 

0. 00344 



44 

7 

3.90 

0.00324 

3.80 

0.00336 

45 


3.76 

0. 00340 



46 


3.99 

0. 00423 



47 

9 

3.80 

0.00442 

3.80 

0.00419 

48 


3.63 

0.00392 



49 


4.00 

0.00536 



50 

15 

3.72 

0.00537 

3.90 

0.00543 

51 


3.98 

0.00558 



52 


3.78 

0.00596 



53 

21 

3.81 

0.00610 

3.79 

0.00615 

54 


3.78 

0.00640 






rable 5. Experimental Results for String Mass Concept with Equidistant 
Lead Shots and with Jam Nuts 


Number of 
lead shots 


Frequency 
f , Hertz 
e 


Dampi ng 
ratio, c 


Aver age 
f , Hertz 
e 

Averane 

C 

3.93 

0.00228 

3.89 

0.00258 

3.98 

0.00462 




























Table 6. Effect of Axial Force on Member Frequency and Damping Ratio 


Test 

No. 

Axial Force, P 
lb. 

Frequency 
f , Hertz 
e 

Dampi ng 
ratio, c 

Average 

f , Hertz 
e 

Average 

c 

64 


6.78 

0.00106 



65 

-155.0 

6.37 

0.00060 

6.64 

0.00083 

66 


6.78 

0.00083 



67 


■RH 

0.00104 



68 

-19.8 

Ilgpa 

0.00103 

4.20 

0.00104 

69 



0.00107 



70 


3.83 

0. 00110 



71 

-4.5 

3.85 

0.00105 

3.83 

0.00105 

72 


3.82 

0.00100 



73 


■PM 

0.00517 



74 

+10.8 


0.00564 

3.42 

0.00526 

75 


mEsm 

0.00499 



































I 


lb 


fable 7. Experimental Results for Polyethylene Tubing Concept 
with Jam NLits 


Test 

Tubing 

Frequency 

Damping 

Average 

Average 

No. 

Arrangement 

f , Hertz 
e 

ratio, c 

f , Hertz 
e 

C 

76 


3.96 

0.00239 



77 

(a) 

3.92 

0.00208 

3.95 

0.00220 

78 

3.96 

0.00213 



79 


4.01 

0.00238 



80 

(b) 

3.74 

0.00195 

3.93 

0.00205 

81 

3.98 

0.00183 



82 


3.72 

0.00199 



83 

(c) 

3.66 

0.00202 

3.68 

0.00207 

84 

3.64 

0.00221 




f 

i 
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Table 9. Experimental Damping Ratios for Chain Danplng Concept with Jam tots 


Test 

No. 



Damping ratio, c 




0 sec 

1 sec 

2 sec 

3 sec 

4 sec 

5 sec 

6 sec 

85 

0.00312 

0.00263 

0.00241 

0.00283 

0.00291 

0.00531 

0.00341 

86 

0.00301 

0.00271 

0.00269 

0.00277 

0.00287 

0.00299 

0.00349 

87 

0.00295 

0.00293 

0.00240 

0.002 70 

0. 00292 

0.00300 

0.00343 

88 

0.00320 

0.00309 

0.00317 

0.00355 

0.00419 

0.00432 

0.00441 

89 

0.00317 

0.00312 

0.00318 

0.00349 

0.00429 

0.00429 

0.00429 

90 

0.00298 

0.00332 

0.00330 

0.00358 

0. 00424 

0.00432 

0. 00456 

91 

0.00391 

0.00419 

0.00409 

0.00423 

0.00479 

0.00539 

0.00613 

92 

0.00401 

0.00423 

0.00408 

0.00418 

0.00492 

0.00542 

0.00529 

93 

0.00416 

0.00421 

0.00416 

0.00416 

0.00475 

0.00539 

0. 00613 

94 

0.00448 

0.00471 

0.00477 

0.00512 

0.00579 

0. 00639 

0. 00310 

95 

0.00471 

0.00461 

0.00481 

0.00503 

0.00583 

0.00621 

0.00801 

96 

0.00464 

0.00498 

0.00509 

0.00512 

0.00584 

0.00665 

0.00792 

97 

0.00539 

0.00567 

0.00601 

0.00659 

0.00701 

0. 00866 

0.01130 

98 

0.00562 

0.00522 

0.00613 

0.00671 

0.00712 

0.00849 

■ 0.01090 

99 

0.00574 

0.00534 

0.00583 

0. 00663 

0.00717 

0.00833 

0.01100 

100 

0. 00633 

0.00659 

0.00741 

0.00912 

0.01170 

0.01731 

0.02330 

101 

0.00639 

0.00658 

0.00753 

0.00S91 

0.01230 

0.01693 

0.02370 

102 

0.00721 

0.00713 

0.00756 

0.00907 

0.01430 

0.01745 

0.02530 

103 

0.00879 

0.00912 

0.01020 

0.01320 

0.02310 

0.02790 

0.04170 

104 

0.00891 

0.00911 

0.01060 

0.01330 

0.02290 

0.02830 

0 .04 110 

10b 

0.00736 

0.00377 

0.01130 

0.01000 

0.02030 

0.02990 

0. 03990 


















Table. 10. Experimental Average Damping Ratios for Chain Damping Concept with Jam Nuts 


Chain 

length, 

feet 



average damping ratio. 

C 



0 sec 

1 sec 

2 sec 

3 sec 

4 sec 

5 sec 

6 sec 

1 

0.00303 

0.00279 

0 00250 

0.00260 

0.00290 

0.00310 

0.00346 

2 

0.00312 

0.00317 

0.00325 

0.00354 

0.00424 

0.00431 

0.00442 

3 

0.00403 

0.00421 

0.00411 

0.00419 

0.00482 

0.00540 

0.00620 

5 

0.00461 

0.00476 

0.00489 

0.00509 

0.00582 

0.00641 

0.00301 

7 

0.00558 

0.00541 

0.00599 

0.00664 

0.00710 

0.00866 

0.01100 

9 

0.00699 

0.00676 

0.00750 

0.00903 

0.01276 

0.01723 

0.02423 

10 

0.00852 

0.00900 

0.01070 

0.01320 

0.02210 

0.02870 

0.04090 
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Table 11. Efficiencies of Various Damping Concepts 


Test 

NOS. 

Details of damping 
devices 

“g. 

Average 

5 

Efficiency 
Index, n 
1b-sec 2 /ft 

1-3 


0 

0.00000 

0.00243 

_ 

4-6 


1 

0.00103 

0.00278 

10.94 

7-9 

Number of lead 

2 

0.00207 

0.00353 

17.11 

10-12 

shots at midheight. 

3 

0.00310 

0.00395 

15.78 

13-15 

without jam nuts 

5 

0.00517 

0.00429 

11.58 

16-18 


7 

0.00723 

0.00466 

9.93 

19-21 


9 

0.00930 

0.00519 

9.55 

22-24 


15 

0.01550 

0.00601 

7.43 

25-27 


21 

0.02170 

0.00693 

6.67 

28-30 


0 

0.00000 

0.00103 

— 

31-33 


1 

0.00103 

0.00145 

13.12 

34-36 

Number of lead 

2 

0.00207 

0.00227 

19.28 

37-39 

shots at midheight. 

3 

0.00310 

0.00263 

16.61 

40-42 

with jam nuts 

5 

0.00517 

0.00281 

11.08 

43-45 


7 

0.00723 

0.00336 

10.37 

46-48 


9 

0.00930 

0.00419 

10.93 

49-51 


15 

0.01550 

0.00543 

9.11 

52-54 


21 

0.02170 

0.00615 

7.59 

55-57 

Number of Equi- 

2 

0.00207 

0.00228 

19.44 

58-60 

distant read shots 

3 

0.00310 

0.00258 

16.10 

61-63 

with jam nuts 

21 

0.02170 

0.00462 

5.31 

76-78 

Polyethylene 

(a) 

0.30480 

0.00220 

0.12 

79-81 

tubing 

(b) 

0.16800 

0.00205 

0.20 

82-84 

concept 

(c) 

0.27300 

0.00207 

0.12 

85-87 

. 

1 

0 • (bOOO 

0.00291 

1.21 

88-90 


2 

0.10000 

0.00372 

0.86 

91-93 


3 

0.15000 

0.00471 

0.79 

94-96 

Chain length. 

5 

0.25000 

0.00565 

0.71 

97-99 

ft. 

7 

0.35000 

0.00719 

0.56 

100-102 


9 

0.45000 

0.01200 

0.79 

103-105 


10 

0.50000 

0.01900 

1.15 
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(b) 


Fioure 1. Lightweight Load Support Vibration Damper (Ref . 7). 
















DEFLECTION. INCHES 


Envelope 


SEMI -ANALYTIC 
FINITE-DIFFERENCE 




S ’ll 


a;\a; 


TIME. SECONDS 

Comparison of Semi -Analytic and Finite-Difference Deflection-time plots 
for c = 0.0060. 
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Figure 11. Spatial Initial Imperfection. 
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Figure 18. A-t plot for Test 1C, Table 3, String-Mass Concept without Jam Nuts. 
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Figure 21. A-t plot for Test 19, Table 3, String-Mass Concept without Jam Nuts. 
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Figure 27. A-t plot for Test 37, Table 4, String-Mass Concept with Jam Nuts 
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43, Table 4, String-Mass Concept with Jam Nuts. 
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Figure 30. A-t plot for Test 46, Table 4, String-Mass Concept w ith Jam Nuts. 



DEFLECTION. INCHES 


.riSV5 5 

ml; is 11 


4 

4 fm 


Y : '\ i J • •:> 

r p .1 . t. . i- 


! -■ : 1 r : 


m 

m 

; r >.] p. 


ail i 


um 


:2 i-Li- 





L; :.T L> ' -V.i s r *£.\ t LI **-; ) H JJ 4;’ U * •• 


TIME. SEC. 

Figure 31. A-t plot for Test 49, Table 4, String-Mass Concept with Jam Nuts. 
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Figure 32. A-t plot for Test 52, Table 4, String-Mass Concept with Jam Nuts. 
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Figure 33. Effect of Connection Friction on c-n Relationship 
with Lead Shots at Midheight. 
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Figure 35. &-t plot for Test 58, Table 5, String-Mass Concept, Equidistant lead Shots 
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Figure 41. A-tj>lot forest 73, Table 6, Effect of Axial Force 
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Figure 42. 


A-t plot for Test 76, Table 7, Polyethylene 
Tubing Concept (a). 
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Figure 43. *-t plot for Test 79, Table 7, Polyethylene 

Tubing Concept (b). 
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Figure 44. 4-t plot for Test 8 

Tubing Concept (c). 
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Figure 45. Details of Chain Links. 
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Figure 46. 


A-t plot for Test 85, Table 8, Chain Damping 
Concept (1 ft. Chain). 
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Figure 55. Comparison of Experimental and Semi-Analytic A-t 
plots for Test 52, String-Mass Concept. 
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Figure 56. 


Comparison of Experimental and Semi-Analytic A-t 
plots for Test 76, Polyethylene Tubing Concept. 
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Figure 57. Comparison of Experimental and Semi -Analytic A-t 
plots for Test 103, Chain Damping Concept. 
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Figure 59. Comparison of Experimental and Finite-Difference 

&-t plots for Test 76, Polyethylene Tubing Concept. 
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Figure 60. 


Comparison of Experimental and Finite-Difference 
A-t plots for Test 103, Chain Damping Concept. 
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Figure 61. Efficiency Index Relationships. 
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APPENOIX A 


; Computer Programs 

j. Computer programs were writtei. based on the semi-analytic and the fin- 
te-dlfference procedures described In Chapter II. The following are the 
nput data for the two computer programs given in this appendix and naned 
PARTA," and "Part B,*' respectively: 

M Number of nodes along the member length = 12, 

2 

E Young's modulus = 30,000 kips/irt, 

i* 

MI Moment of inertia = 0.00213798 in, 

L Length of the member * 144.0 in, 

K Rotational stiffness = 2.230 kip-in/rad, 

DO Initial member deflection = 0.366 in, 

2 2 

Rho Mass per unit length = 651.0 x 10 -10 kip-sec /in , 

Zeta Damping Ratio. 

Actually the M value is not needed in Part A but is included keep 
the input data for both Parts A and B identical. In Part B, the time Incre- 
ment At can be adjusted as needed by modifying the statement numbered 121 
(T=T+0.001) in which At is taken as 0.001 secs. For example, T = T + 

0.0005 may also be used. Also, to obtain the A-t output for t > 0.5 
secs., the fifth statement above the statement numbered 121 may be modified 

(from DO 7 NN = 0,500,1 for t = 0.5 secs., to say, DO 7 NN = 0,7000,1 for 
t = 7 secs). 

The output for each of the programs is: 

Damping ratio 
Rotational stiffness 

Time and deflection at every 0.005 seconds from the start. 
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PART A 

S£HI_ANALYTICaL COMPUTER PROCECURE 

xxxxxxxxxx 7 x:xxxsiix::x::::x:xxxxx:x:xs;xtxxx::x:: 

• ♦ PRCJECTX DYNAMICS OF A P APT I ALLY PESTRUNGP CtiLD. M N'** 
MMMM*** PROF X 0K.2IA RA2ZAO * •*•*•♦•* * 

**'«***«*** 3 Y S RAJEN0R4 KUMAR EKHELlKAP. 


DEFINING OF YARlAdLcS 


DOUBLE PRECISION A (20, 400 » , X ( 20 , 1) , Y (20 > * ^F («* CCC ) » j S ( l?9J 
DOUBLE PRECISION CC( 5001. '*J<20 ) »HK( 2tl , CO A 4 »L AMC, P C R * H 
DOUBLE PRECISION B l, !>2.33 .P4.55 » c6, 37 ,23 ♦ a9 , £ » • I . PC * 
DOUBLE PRECISION L»DANl,DAM2»Pu»Cl,0*,$»2.W‘'»')C»P 
OCUBlE PRECISION ZZZ.YYY.PtE.ZETA.K.CC.CG 
OQUbLE PRECISION HG( 4COO J »R 1. RZ , *3. LLU t G«i 
INTEGER M.I ,.!,KK,NN,N, IKK 


READING OF VARIABLES 


DATA M.E.MI »L»K»U0/12»3G0G0.»0.0C21.47 < »e»i4*,»2.2j,. 3o6/ 
DATA ZETA.R3W/0.5000,fc51. CC37E-10/ 


CALCULATIONS OF DAMPING CCEFF 


* 

♦ 

# 


PIE -3.1415927 
MRITF. (22»12 3)ZETA 

FORMAT! 0 AMP ING COEFFICIENT IS‘Fb.3) 

WRITE (22» 1251K 

F ORMAT C • THE ROTATIONAL STIFFNESS IS*F8.3> 

WRITE (22.1231 

FORMAT! • THE TIME ANO THE CORRESPCND INC OEFLECT IONS AkE*> 
OAMl-(( 12 .00* ( (PIE*G*MI)**2 )* ! P I E ** 2 ) ) / L* *2 ) 

♦ !32.D0*E*rtI*K*L* t (5.00* (P IE**2 ) ) I / (2.DC*L**2 ) ) ♦ 

(3.0 0M(K*LI**2I*4.D0*(PIE**2>/CL**2»> 

0AM2- ( 12. DO* { (PIc*E*MI ) ** 2 ) 1*1 32.00*t*rtl*K*L) 

♦ C (9 .DC* (K*L1**2I /4 .00) 


LAMC- C < PlE/L J **2 )• (DAM1/0AM2 J 
OM- (£*M1»LAM0/KGW) 

CDAK-2.00*Z6TA*QSQRT(E*MI ♦RCW^LAHO) 
1-0 
T-0. 


ORIGINAL PACE 13 
OF. POOR QUALITY 


00 14 1-1,100.) 
ICOUNT-ICGUNT*! 

IFt ICOUNT .GT.2500) I COUNT* 2 5 CO 


T-T+.C05 


Kl-OSCRTt 4.()0*0M-{CDAM**2/Rnvi**2 ) ) 

R3-1.0OM (K*L)/(2.00*G*PIE*MI) I 
R 2-00 /R 3 

WG( I)“DO* (DExP{-(COAM*T/( 2.no*K(l«) )) ) ♦ < OCOS ( R l* T/2 . DO > 
♦ ( X CD AM/( ROW* R1 )) ♦OSIN (R1 *1/2.00 )) ) 

WRITc(22,*)T,WGCH 

CONTINUE 

STOP 
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DAMPING RATIO (ZETA) LS 0.500 
ROTATIONAL STIFFNESS IS 2.230 
TIME AND THE CORRESPONDING OEFLECTI QNS 


ARE 


5.0C0000E-03 
1.000000E-02 
1.5C0C00E-02 
2.0C0000E-02 
2.590000E-02 
3.GCOOOOE-02 
3.5COCOOE-02 
A.CC0000E-02 
4.5C0000E-02 
5.0C0Q00E-02 
5.500 000E-02 
6.0000006-02 
6 .5 COCQOG-O 2 
7.GCQ000E— 02 
7.500000E-02 
8.0C0G00E-02 
0.5COOOOE-O2 
9.0000006-01 
9.500000E-02 


O.ICOOOOO 
0.1050000 
0.1100000 
0.1150000 
0.1200000 
0.1250000 
0.1300000 
0.1350000 
0.1400000 
0.1450000 
0.1500000 
0.1550000 
.0.1600000 
0. 1650000 
0.1700000 
0.1750000 
0.1600000 
0.1050000 
0 . 1900000 
0 • 1950000 
0.2C00000 
0 . 2C5C000 
0.4100000 
0.2150000 
0.2200000 
0.2250000 
0.2300000 
0.2350000 
G. 2400000 
0.2 460000 
0.2500000 
0.2550000 
0.2600000 
0.2650000 
0.2700000 
0.2760000 
0.26C0000 


0.36359234245766419 
0.35675841937089732 
0.34607516684695260 
0.33210549126955565 
0.31539255379189680 
0.29645494226637371 
0.27578269234247022 
0.25383411381816543 
0.23103333186338223 
0.20776883146797169 
0.18439190347081269 
0.16121668308238770 
0.13851997143343834 
0.11654183037001003 
9. 54365 4 1208 840 3 41E -02 
7.552 39193712530 97 E-02 
5. 6790928634 50 8939E-G 2 
3. 9 39 35410853372 97 E-02 
2. 340879164290 37 33E-0 2 
0869791536118727E-03 
-4. 146030531492 2873E-0 3 
—1.568 64 4C478508267E-02 
-2 .5 75 01 34 83 548 33 5 3 E-02 
-3. 43703006781961 44E-0 2 
-4 .159 50 8915775 97 19 E-02 
-4. 748534277012 8493E-0 2 
-5.211240555431 16 61 E-02 
-5. 5 55 604643 85324 34 E-02 
-5. 79024982493764 2 3 E-02 
-5. 924 2631 14752 22 d6E-0 2 
-5. 967026790715 74 5b E-02 
-5. 928C647 11792 9573 E-02 
-5. 8 169037 80 813 74 75 E-02 
-5.6429506647467103E-02 
-5. 4 15 383679065 6145 E-02 
-5 .143059557020 32 45 E-02 
-4. 8 34 43 4663 363 51 25 E-02 
-4 .497500074 3 14 43 37E-0 2 
-4.13972983C4452952E-02 
-3. 768 04 1575 73595 25E-0 2 
-3. 3 88768723073 40 03E-0 2 
-3.00764323265511 lOE-O 2 
-2. t>29 78 30 53 669 49 14 E-0 2 
-2. 25971 6259 700 33 31E-0 2 
-1.90134990379811226-02 
-1.5530 15624 922 296 l E-02 
-1.23243606 1. 53063 87E-0 2 
-9. 2699b 152 1416 39 o4£-0 3 
-6. 4 32 754453 894 37 4 l E-0 3 
-3 .625015351935 4139 e-0 3 
-1. 4 5736 272 12 19 92 90 E-0 3 
6. 6 837 462999 505 17 5E-04 
2. 550764 3163241710 E-0 3 
4. 192384 59707549296-03 
5. 59 36 36 1 422 641 91 2E-0 3 
6. 7773606891 7 72 66 7E-03 
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0.2650000 
0*2000000 
0.2050000 
0.3000000 
0. 3C50000 
0. 3100000 
0.3150000 
0.3200000 
0.3250000 
0.3300000 
0.3350000 
0.3400000 
0 .3450000 
0.3500000 
0.3550000 
0.5600000 
0.3650000 
0.3700000 
0.3750000 
o.3eooooo 
0.3650000 
0.3900000 
0.3550000 
O.4COG0OO 
0.4C50000 
0.4 100000 
0.4150000 
0.4200000 
0.4250000 
0.4300000 
0.4350000 
0.4400000 
0.4450000 
0.4500000 
0.4550000 
0.4600000 
0.4650000 
0.4700000 
0.4750000 
0.4800000 
0.4650000 
0.4900000 
0.4950000 
0.5000000 


7. 7 384 70 30 46350 44 4E-0 3 
8. 4935903668663 86 7E-03 
9. 05 5722 7970 485 30 0E-03 
5. 438921 51 44 204 26 IE-03 
9. 65 79995606608 26 7E-03 
9. 72 025104 4946584 IE-03 
9. 6652050022 53561 6E-03 
4843914370596397E-03 
9.2011 462424 496 712E— 03 
e. 8304408429 123 143E-03 
e. 38671090671 31 52 7E-03 
7.8837478106821786E-03 
7.33458271550d4403E-03 
6. 75 1404 12 185 1706 7E-03 
6. 1454936251 132 05 4 E-03 
5. 5271004676211 67 8E-03 
4. 9056 13 3991 329 06 7E-03 
4.289748297bl35574E-03 
3. 60 63 499698 399 51 7E-03 
3. 1020065427657 102E-O3 
2.5421548b98321568E-03 
2. 0113 154079338930E-03 
1. 513135069309853 4E-03 
1.05043661474398886-03 
6. 25 27 322838021 26 7E-04 
2. J898699772005989E-04 
-1 .0772988708 17 09 53 6— 04 
-4. 147723046560 15 43E-0 4 
-6 .025 58735 29 02 1C 68E— 0 4 
-9.119680447' • 5 0735E-O4 
-1.104277307 30 2 69 43 E-03 
-1. 2611013551786858 E— 03 
-1. 384334644425 03 37E-0 3 
-1.476095931 382 83 03E-0 3 
-1.5 3 8676163 2 50 70 15 E-03 
-1.57446 99 0051430246-03 
-1.5 86030508044 46 25 E-03 
-1.575 829277000 35 09E-0 3 
-1.546418571028 1397E-0 3 
-1.5 CO 29902 7 9 09 26 10 E-0 3 
-1.4 39910 792 667 3906E-0 3 
-1.367600706594 87 87 E-03 
-1. 28564133822868986-03 
-1.19613370004031216-03 


CCC PART e 

CCC FINITE DIFFERENCE COMPUTER PROGRAM 

CCC ♦♦ PROJECT 3 DYNAMICS OF A PARTIALLY PESTRIANED COLUMN *♦ 

CCC ***«♦♦*♦♦♦♦ prof : OR* 21 A RAZZ AQ ♦*♦♦♦•»♦**♦*♦♦*♦•♦**♦♦* 

CCC *♦**♦*♦#♦♦♦ BY * RAJENORA KUMAR EKHELlKAR ♦** ♦♦♦****♦♦*• 

CCC *•* «* ♦♦♦*♦**♦*•♦*♦♦♦#♦♦**♦****•*•♦<>*♦»*****♦♦ ♦♦♦**♦♦***• 


CCC 

CCC 

CCC DEFINING OF VARIABLES 

CCC 

CCC 

OCU8LE PRECISION A (20, 400 J » X( 20. I ) , Y 1 20 ) » JF ( 4 COO) , SS ( 150) 
OCUBLE PRECISION CC< 500J . W J (20 » »NK( 20 ) .CDAM.L AMD.PCRtH 
OOUBLE PRECISION 8 It B2 .83 » B A . B 5. B6» B 7.88, 39 »E ,MI ,K0 A 
OQUBLE PRECISION L *0 AMl.D AM2 ,P0 , Cl, 02 ,S ,Z ,WR, NO, P 
OCUBLE PRECISION ZZZ, YYY. PIE.ZETA.K.UC.CQ 
OOUBLE PRECISION WG( 4000) .R1.R2.R3.LGU.QM 
INTEGER M,i,J,KK,NN,N,IKK 

CCC 

CCC 

CCC REAOING OF VARIABLES 

OCC 

CCC 

DATA N,E»MI,L,K.0D/12,300C0.,0.0G214793»144.,2.23».366/ 
OAT A ZETA.R0H/0.5000* 651. 00876-10/ 

CCC 

CCC 

CCC CALCULATIONS OF DAMPING CGEFF 

CCC 

CCC 

PIE -3.1415927 

DAM1«({12.D0*((PIE*£*MI)**2)*(PIE**2))/L**2 ) 

♦ ♦ (32.03*E*MI*K*LM(5.00MPlE**2)))/(2.00*L**2) )♦ 

♦ I3.00*(IK*L)**2)*4.00*(PIE**2)/(L**2)) 
DAM2-(12.00*((PIt:*E*Ht)**2))M32.D0*E*Ml*K + L) 

♦ ♦( <9.00*<K*L)**2)/4.00) 

LAMG-((PIE/L)**2)*(0AM1/DAM2) 

OH- (E*M I* LA MO /ROW) 

T-0.0 
DT-C.001 
B3-R0W/ (DT'**2 ) 

H-L/12.00 

IKK-0 

Bl-£*MI/( (H)**4) 

B5-((E*MI/(H> )-K/2.00)/( (E*MI/(H) >*K/2.00» 

P-0. 

INN-BOO-1 
ICOUNT-O 
DO 7 NN-0.500.1 
IFINN.NE.QJGO TO 121 
T-0.0 
GO TO 10 
T YP E* »NN, T 
121 T-T*0.001 

IF (T .EQ. C. 00 1 ) GO TO 20 
CDAM-2.00*ZETA*0SQKT (E*MI ♦ROW*LAMO> 

U4-CDAH/DT 
Cl— 1. 00/(03*04) 





ccc 
ccc 
ccc 
' ccc 
ccc 


31 


30 


40 


50 

CCC 

CCC 

CCC 

CCC 

CCC 

10 


70 


CCC 

ccc 

ccc 

ccc 

ccc 

20 


75 


CALCULATIONS OF CONSTANTS UJ i*B2 »B3. •• IFOR MATRIX 


B2**P/ (H**21 
DO 31 IK-1.H-1*1 
ACIK»IK)-0.0 

A ClU»-6- 00*81-2.00*02-2. 00*B3-Bl*b5-B4 
ACM-1»M-1)-A(1,1) 

DQ 30 IK-2.M-2.1 _ 

A l IK* IK)-6.DO*Bl-2.DO*32-2. 00*B3-B4 
CONTINUE 

DO 40 I K» 1*H— 2 ♦ l 
A ( I K* IK*i)-B2-4. 00*81 
A(IK*i»IK)-A( IK »IK*1 ) 

CONTINUE 

00 50 IK-l*M-3»l 

AC IK*2*IK)»B1 

AC IK* IK*2 )-AIIK*2» IK J 

CONTINUE 

IFIT.GT.DTJGO TO 90 


CALCULATIONS OF INITIAL DEFLECTION! T-O.O) 


DO 70 I“1*M— 1* 1 
B8»CK*L)/C4.00*PIE*E*NI) 

YYY-FLOATCI) 

Z22-FLOATCM) 

7 . p xe#yyy/zzz 

WK( I )»COD/C 1.00 *96 *2 .00) ) ♦ (OSIN ( Z )*B6* t l.DO-DCOS ( 2*Z ) ) 1 

CONTINUE 

PRINT** WK (6) 

GO IQ 7 


CALCULATION OF FIRST TIME OEFLECTIOMT^O.OCl) 


B5-C(E*MI/12.D0)-K/2.00)/(CE*MI/12.00)*K/2.0O) 

B6— l (DT**2)*E*MI)/< 2.00* ROW* ( 1 12.0 01**4)) 
b7 «_(P*(DT**2I ) / ( 2 .DO* ROW * ( 12. 00)** 2 I 
B8«-(P*(DT**2)> /C 2.0 0* ROW ) 

DO 75 1*1. M— 1* 1 

YYY-FLOAT (I) 

ZZZ-FLQATtn) 

Z-P IE*YYY/ZZZ 

Y C I >«-DD*UPIE**2»*SIN(Z) )/(L**2) 

CONTINUE 
J-l 
I "1 

W J( I J-(6.DO*Bfa-2.DO*B7*l.DO)*HKI I )-4 .00 ♦ B fa * «K ( 2) 

* ♦ (B7*WK(2 )*3b*hK{3 J + U8*Y( I ) -B5 * Bb * *K Cl) > 

I “2 

WJl I) = t 6.00*Bfa-2. DO* B7*l. DO )*WK( 2 I*C 137-4.00*86 1 *HK(ll 

♦ ♦(B7-4.D0*&6)*MK(3)«-Bb*WK{4)'*-8B*Y(l > 

WJCiD-WJCll 
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I 


WJ(10)-WJ(2) 

00 80 I»3tM-3tl 

WJCI)"B6*WKlI--2)*{B7-A.00*B6)*wK(I-l)'M6.D0*B6-2.D0*B7 
♦ l.CO)*HK(IJ*( B7-A >00*86) ♦ WKI I * 1 ) ♦() 6* WK ( I *2 ) ♦ B8*Y ( I ) 
CONTINUE 
J-0 

UR I TE (21» 111) K 

F CRPATI • THE ROTATIONAL STIFFNESS IS»,F8.3> 
URITEI2I.1I212ETA 

FQRPATC * THE DAMPING COEFFICIENT IS SF8.3) 

WRITE (21, 113) 

FURMa'TCTHE TIME AND THE CORROSPOND ING DEFLECTIONS ARE') 
GO TO 7 


THE DEFLECTIONS AT TIME GREATER THAN THE FIRST INCREMENT 


IKK* IKK* I 
DO 100 I-1«M-1,1 
YYY-FLOAT(I) 

Y(I1— 0D*C(PIE**2)*DSIN(P IE ♦YYY/12. DC) 1/20736.00 

SUM-0.0 

DO 122 KK-1 til 

SUM»SUM*A(ItKK)*WJ(KK) 

CONTINUE 

HFCI)-C1*SUM»D2 *WK (I)+P*Cl*Y(I) 

CONTINUE 

DO A3 A 1-1* M— I* 1 
NKC I) -W JC I) 

WJ(I)-HF(I) 

CONTINUE 

IFt Ii'K.EQ.5)WRITE(21,*)T,WF(6) 

I F l IKK.EQ.5IIKK-IKK-5 

CONTINUE 

STOP 

END 









6 

1 

1 

2 

2 

3 

3 

4 

4 

5 
b 

6 
6 
7 

7 
a 

8 
9 
9 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


ROTATIONAL STIFFNESS IS 2.23 

damping ratio (zeta) is 0.5 

TIME AND THE CORROSPQND INC DEF 


0C0000E-Q3 
100000E-02 
600000E-02 
1COOOOE-02 
600000E-02 
lOOOOOE— 02 
6C0000E-02 
1C0000E-02 
6G0000E-0Z 
100000E-02 
bOOOOOE-OZ 
1G0000E-02 
bCGOOGE-02 


0.36054476463108740 
0.34857072755345545 
0.33162699131 269795 
0.31530570391608760 
0.30050586557196719 
0.28609830892547833 
0.27068924323488135 
0.25393064779027415 
0.23405013350048578 
0.21162899446399367 
0.18647207550126029 
0.16164113546252522 
0.13868620961410036 


• 1000'JIE 
.60000 IE 
• 1C0001E 

• 6 0000 IE 
•1Q0001E 
• 6 C 0001 E 
.1010000 
.1060000 
.1110000 
.1160000 
.1210000 
.1260000 
.1310000 
.1360000 
.1410000 
.1460000 
.1510000 
•15o0000 
.1610000 
.1660000 
.1710000 
.1760000 
.1810000 


02i 

02 , 

02 , 

02 , 

02 , 

02 , 

2 . 

5 . 

-7 

-1 

-2 

-3 

-3 

-4 

-4 

-4 

-5 

-5 

-5 

-5 

-5 

-5 

-4 


0.11752472410034994 
9.9 00 79 31338 57 68 75 E 
8.3008°74406320450E 
6.7 72 0775151456077 E 
5. 1912716058 67 4933E 
3.5830642709889558E 
04 20828203 13928 7E-02 
59 02 220 3 62 88 907 8 E-03 
.6C91531849992118E-0 
.82676 1919632 9048E-0 
•6108876511530570E-0 
.227989077 325 2890E-0 

• 74825393315633 84E-0 
•178 43 3541 668 9278 E-0 
.5 86510824086 2084E-0 
. 984 4546924 17 00 27E-0 
. 338975401961 4770E-0 
.5 5339809 1524 3934 E-0 

• b2842314?9 307660E— 0 
.584070062 3 6698 54E-0 
•427204511954 56 39 E-0 
. 188 340006 7909164E-0 
•9 3 199 3 09 1030 04 11 E-0 


.1660000 
.1910000 
.1960000 
.2010000 
.2060000 
.2110000 
. 21 o 0000 
.2210000 
.2260000 
.2310000 
.2360000 
.2410000 
.2460000 
.2510000 
. 2 5bC000 
.2610000 
.2660000 


-4.7C8335080527 75 5 6 E-0 
-4. 491 17585 4 698 48 64E-0 
-4. 263 6078 26607 49 88E-0 
-4 .012612471727 86 13 E-0 
-3. 7 32 102554610 3294 E-0 
-3. 390 404439748 88 132-0 
-3.0 0479 34 97 2 12 74 46E-0 
-2. 6 19 22 54 52 3 73 50 46E-0 
-2. 2 56 33452 2 4d8 45456—0 
-1.9 25 106469 463 4b59E-0 
-1.6262909 56664 27 53E-0 
-1.37268 79 5 76 35 91 36 E-0 
-1.1 34 9695 8 3b 27 8 2 38 E-0 
-8.9 17 Id 10 3 8 827 8820 E-0 
-6. 4 16 CO 64 99 2 86 05 01 E-0 
-3.973774022032 12 39 E-0 
- 1.651 12976 4 4 61 76 66 c-0 


.2710000 

.276C000 

.2810000 


4.967323527O985887E-04 
2. 2552525356390655E-03 
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. 5970 635743793 30 7E -03 
. 45 573 8 72 12 990 89 8 E-03 
. 1707471795620414E-03 
.8213707394299620E-03 
.46455255381 75 598E -03 
,055356076«2b7363E-03 
.4686949717737365E-03 
. 63 90 865653 595 87 2E -03 
.6312933065495976E-03 
.4358 8761 12904 440E-03 
, 1048616373026998E-03 
, 7148 725016753 97 IE-03 
,3655721062 142660E-03 
, 0463 2351 1302 4 41 8E-0 3 
, 703511 342 7 9 10 72 4E-0 3 
.3352377532903728E-03 
>9164013869270539 E-03 
, 419085 8062 857680E-03 
, 833721799720807 3E-03 
.2344 126759690 29 7E-03 
. 665757 8003 13061 5E-03 
. 1426644653 173 22 5E-0 3 
, 6689733915 879 129E-03 
.2597661755 80908 5E-03 
, 892174 8153 578975E-03 
, 5 159 370075395 46 IE-03 
. 1294633324869032E-03 
. 4361 579156653 096E-04 
.7854482535022 735E-04 
, 46 78715077558 18 7E-05 
! • 567659122 3 2038 29E—0 4 
». 7882 1109142920856-04 
>•4 849407938908323c— 04 
r. 869012876013 9707E-0 4 
>•073816178981 99 87 E-0 4 
L. 011086105 397 50 71 E-03 
L.1143984379016757E-03 
L.211 97 059120673 b4 E— 0 3 
L.28o6509999497057E-03 
L.32422342797723b7E-03 
L. 3 30 38 6004 474 4655 E-03 
L. 3095775492935182 E-03 
L. 2637225988351805 E-03 
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► Instrumentation and Test Setup 

V 

) 

; A brief description of the Instrumentation and test setup used for 
inducting the vibration experiments described In Chapter III Is given here, 
gure 10 shows a system consisting of a probe, vibration Instrumentation, 
id a deflection-time curve plotter schematically. A brief description of 
■ ie system components is described In Section B.l of this appendix, while 
he procedure for setting up the test apparatus Is outlined in Section B.2. 

.1 Instrumentation 

The principal components of the instrumentation and their respective 
unctions are described below: 

1. Probe or sensor, and oscillator demodulator. 

This component of the system Is manufactured by Kaman Measur- 
ing Instruments (Model KO 2300-10C). The probe consists of a vari- 
able impedence bridge with an active and a reference coil. The 
variation of Impedence results from the eddy currents Induced In 
nearly conductive surfaces. The probe gives exceptional accuracy 
for non-contacting measurements of conductive surface motions. It 
has a high capability for linear calibration within ± 1.0 percent 
deviation with almost any metal or alloy, whether flat or curved. 
Adjustment controls are provided such that the output voltage is 
readily calibrated so that It Is directly proportional to the 
displacement and has a wide frequency response (0.20 KHz, ldb 
point; and 0.5 KHz, -3dB point). 

The stability Is affected by the target material such as a 
magnetic material. It works very well when the target Is near the 
face of the sensor. The movement of the cable connecting the 
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sensor and the oscillator demodulator, the dielectric constant, 
sensitivity to noise, temperature, hunidity, and the magnetic field 
are all greatest when the target Is at Its full scale displacement. 
Linearity Is not observed when the distance of the target from the 
sensor is within 10 percent of the full scale displacement. The 
probe converts the displacement of the conductive surface into an 
output voltage. 

Galvanometer amplifier and data graph. 

The main function of the galvanometer amplifier is to accept 
the signal from the transducers, that Is, from the probe and the 
oscillator demodulator which translate mechanical variation into 
electrical impulses, and Is conditioned before sending the signal 
Into the data graph. The amplifier is uniquely designed to facili- 
tate the recording of the signal voltage in a light beam oscillo- 
graph. The amplitude of the output signal voltage that can be 
accommodated by the amplifier ranges from 0.1 to 300 volts with 0.5 
to 3 inches peak galvanometer deflection with proper setting of 
volts/in., gain, and the offset knobs. 

The oscillator demodulator removes the carrier frequency from 
the envelope and extracts only the signal frequency that is to be 
fed into the galvanometer amplifier. 

The data graph consists of five knobs which control the output 
graph. 

a. Power on/off, to be connected to a 230V, 60 Hz power 
supply. 

b. Light intensity knob, to control the intensity of the 
light from the recording lamp to the galvanometer. 
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c. 



Paper speed knob, to control the speed of the photographic 
film. The speed range is 0.1 to 100 in/sec., selected in 


| steps of 0.1, 0.2, 0.5, 1.0, 2.0, 5.0, 10, 20, 50, and 

I 100 In/sec. 

d. Gridline knob, consisting of two push button switches to 
| provide "fine" (0.1 In.) and "coarse" (0.5 in.) gridline 

spaclngs, respectively. 

f 

I e. Transport knob, consisting of two Interlocked push buttons 

! to actuate the paper drive. One Is a locking type used to 

activate a continuous run. The other provides a "stop" 
function as well as a momentary run or "jog" function. 
Electrical impulses produced by the transducers may be fed into an 
amplifier where the signal is conditioned before being Introduced into the 
oscillograph. The signal Is then applied to a highly sensitive galvanometer 
in the oscillograph. Light from a mercury arc or xenon lamp is reflected by 
? the galvanometer mirror through a simple optical system which produces a 
spot of light on a moving strip of photographic material. By using various 
; combinations of light intensity and record speed, clear legible records are 
‘ attainable over a wide range of apnli cat Ions. 

• B.2 Test Setup 

; Tiie steps given below are followed for setting up the test apparatus: 

; 1. Connect the instruments as shown In the schematic circuit diagram 

as shown In Fig. 62. 

2. A 120V, 60Hz power is supplied to the galvanometer amplifier, 
datagraph and the oscillator demodulator. 

3. Pull the specimen by approximately 0.75 in. from the p r cbe and 
adjust che three knobs of the galvanometer amplifier such that the 
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maximum deflection on the datagraph is approximately one-and-a-half 
Inches from the center. Adjust the volts/in and the gain knobs to 
their maximum values, and the offset knobs at zero to obtain the 
desired output. 

4. Adjust the datagraph knobs by setting the paper speed to 2 in/sec 
for low frequencies, and select the light intensity to obtain an 
appropriate thickness of the line describing the A-t plot. Fine 
gridlines are usually preferred for a good quick reference. 

5. The specimen is pulled away or towards the probe by means of a 
nylon string. The specimen is then released by burning the string 
with a matchstick to Induce natural vibration. Before burning the 
string, the datagraph switch button "run" is pushed to start. 

6. The datagraph is made to run for 8-10 seconds or 20 inches of the 
plotting paper to get the desired output. 
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